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Abstract. We consider a non-local isoperimetric problem with a repulsive Coulombic term. 
, In dimension three this corresponds to the Gamow's famous liquid drop model. We show that 

whenever the mass is small the ball is the unique minimizcr of the problem. The proof is based 
f***^ ■ on a strong version of a quantitative isoperimetric inequality introduced in [4] 

(N 
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1. Introduction 

m 

We study a ground state problem with a volume constraint 
OH" (1.1) min(p(W f f^0 dydx: \F\ = m 

•<r ■ \ jR n Jr" \ x — y\ 
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where F C K™ is a set of finite perimeter and xf denotes its characteristic function. Here P(F) is 
the perimeter of F and the term with the Newton potential is called non-local part of the energy. 
It acts as a repulsive force and prefers sets to be disconnected. 

The most important case is n = 3, which corresponds to Gamow's famous liquid drop model 
[5]. It was introduced to model the stability of the atomic nucleus and the atomic fission. The 
nature of the problem (1.1) is very fundamental and it appears in many physical phenomena. For 
^ ' more about the physical background and further references see [8] , where the general dimensional 

case is considered. 

In the absence of the repulsive non-local part of the energy, the above problem is the isoperi- 
metric problem. It is well known that then the unique minimizer, up to a translation, is the ball, 
and it is therefore plausible that the ball remains the minimizer of (1.1) when the repulsive part of 
the energy is small. However, due to the different behaviour of the two competing energy terms, 
| the problem is very delicate and even the existence of a minimizer is highly non-trivial. 

The goal is to prove the following stability result for small masses. 

Theorem 1.1. Suppose that n > 3. There exists m n > 0, depending only on the dimension n, 
such that for every m < m n the ball B Tm is, up to a translation, the unique minimizer of (1.1). 

The above result was proven in a very recent paper by Kniipfer and Muratov [7] in dimensions 
n < 7. There the non-local part of the energy is allowed to have a general Riesz potential structure 
Jr- Jk« XF \x-yp V) d V dx with a<n-l, which covers the case of the Newton potential. They also 
considered the problem of unstability for the large masses. A somewhat similar approach to a 
slightly different problem can be found in [3]. 

We stress that this paper was done independently from [7]. Our approach is quite different 
and it is based on a new quantitative isoperimetric inequality introduced in [4]. In fact, the most 
important case n = 3 is nothing but a direct corollary of Theorem 2.1 stated below. We need 
some extra work in the higher dimensional cases but the arguments remain standard. We choose 
to deal only with the Newton potential, to avoid the simplicity of our idea getting buried under 
technicalities. Moreover, Theorem 1.1 is, to best of the author's knowledge, the only result in 
higher dimensions. 

Finally we note that the planar case was studied extensively in [6] with a general Riesz potential 
in the non-local part of the energy. We remark that Theorem 1.1 is not true in the plane, i.e., 
when the non-local part of the energy scales as the Newton potential. This is due to the fact that 
in the plane the Newton potential behaves as log(l/V). Since log(l/r) — > — oo as r — > oo one may 
consider a family of sets consisting on two disks which are further and further from each other. It 
is then obvious, that the value of the problem corresponding to (1.1), would always be — oo. 
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2. Proofs 

By scaling we may write the problem as 
(2.1) min(V) + A m / f f^^ dydx: 

V Jr-Jr" \x-y\ n 



\E\ = IB, 



where A m = (pJTj) an d to is the mass in the original problem. In future we prefer the sets to 

have the measure of the unit ball and therefore we consider the problem (1.1) in scaled form (2.1). 

As mentioned in the introduction, the proof is based on a new type of quantitative isoperimetric 
inequality. Suppose that we have a set of finite perimeter E such that \E\ = \B\\. We measure 
the oscillation of the boundary of E by the assymmctry 

(3(E) := imnjQ^ ^(^-^^(^(a;))! 2 ^"- 1 ^)) ' 

where VB 1 {y)(j^y(x)) = [fEfi an d d*E is the reduced boundary of E. For precise definitions and 
properties of the sets of finite perimeter we refer to [1]. The quantitative isoperimetric inequality 
in strong form reads as follows. 

Theorem 2.1. Suppose n > 2. There is a dimensional constant c n such that for every set of 
finite perimeter E C K n with \E\ = \B\\ it holds 

P(E)-P(B 1 ) >c n (3(E) 2 . 

Notice that by the divergence theorem we may write 

\ [ \v E {x)-v Bl{v) {-K y {x))\ 2 dU n -\x) = \ f \v E ( x )-fll-\ 2 dW-\x) 

1 Jd*E 1 Jd'E \ x ~ V\ 

a*E\ \x-y\J 

= P(E)- f ^±- dU n ~\x) 
Je \x-V\ 

= p(e)-p(b 1 )+ [ I!—l-dn n - 1 (x)-f ^^-dTe-^x). 

Jb 1 (v) \x - y\ Je \x - y\ 

Therefore by defining an assymmctry 

(2.2) 7 (B) := min I f ^— \dx - j ^— \dx\ , 



Bi(y) 



Li- 



Theorem 2.1 can be written as 

(2.3) P(E)-P(B 1 ) >C nl {E), 



for every set of finite perimeter E C l n with \E\ = \Bi\, where C n = . 

We are ready to prove Theorem 1.1. The proof in dimension three turns out to be the easiest. 
This is simply due to the fact that the asymmetry (2.2) scales exactly as the Newton potential. 
The proof also gives a nice bound for the critical mass 7713 > ^ where C3 is the constant from 
(2.3). Unfortunately the proof of (2.3) in [4] doesn't give any explicit bound for the constant C3. 
The higher dimensional cases are proven simlarly by using Lemma 2.2 iteratively. 

Proof of the Theorem 1.1 in the case n — 3. We begin the proof in a general dimension n > 
3 and specify to the case n = 3 only at the end. 

Consider the problem in the scaled form (2.1). Suppose E C R™ is such that \E\ = \Bi\. Since 
the problem is translation invariant we may assume E is centered at the origin, i.e., 

f n ~ 1 j f n ~ 1 1 
l( E ) = / -r-i-dx - / ——dx. 

Jb 1 M Je M 
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Write NL(E) := / RB / Rn x f x { ^Z { 2 } dydx. We need to show that 

NL(Si) - NL(E) < C (P(E) - P{B{)) 
for some dimensional constant C > 0. We organize the terms 
NL(Bi) -NL(E) 

XB 1 (x)x Bl {y) dydx _[ f xe{x) X e(s) dydx 



\x-y\ n 2 ' 7 R «7 R ™ \x-y\ n 

_„ f f XBMiXBi - XE)(X) 

\x — y\ n 
(Xb 1 - Xe)(v)(xb 1 - Xe){x) 



\x — y\ n 2 



dx dy. 



Define 



which solves 



-Av = c n (xbi - Xe) 
for some dimensional constant c„ > 0. Therefore 

(Xb 1 - Xe)(v)(xBi ~ Xe)(x) 



dx dy 

\x — y\ n 1 
= -— [ v(y)Av(y)dy 

C n JR3 

= J- f \Vv\ 2 dy > 0. 

c n JR3 

This and (2.4) yield 

(2.5) NL(B 1 )-NL(E)<2 / / XB 1 (y)(XB 1 -Xe)(x) dxdy ^ 2 f v{y) dy 

Jw Jr." \ x — y\ J Bi 

Assume now that n = 3. Then v(y) — J R „ ^i^rz|p-^ dx. Since v solves — Av = c 3 (xbi — Xe), 
it is superharmonic in B\ and therefore, by the mean value property, we conclude that 

v(y)dy <2v(0) = 2 f ixBl ~ Xe){x) dx = 7 (£). 
Hence, by (2.5) and (2.3), we have 

NL(Bi) - NL(E') = 7(^0 < ^ (P(E) - P(Bt)) 

^3 

which concludes the proof in the case n = 3. 

□ 

The above proof contains all the relevant arguments for the higher dimensional case. Only the 
last step, where we used the mean value property, doesn't generalize since the inequality of the 
type 

v(0) = [ {XB \~ XE 7 ){X) < P(E) - P{B X ) 



is not true when n > 4. Here the notation "<" means that the above inequality is true up to 
a multiplication with a positive dimensional constant. This can be easily seen by considering a 
family of annuli E £ = Br e \ B e , where R £ = \/l + e n . We need to execute slightly more careful 
analysis. To that end we recall some well known definitions and results. 

Suppose that / : M™ — s- R is a locally integrable function. The Riesz potential I a f of /, for 
a € (0, n), is defined by 

(2.6) /q/( , ):=(7 „ q ^_Z|L_ C 
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where the normalization constant is given by 



1 n/ 2iw r(q/2) 

— 7T 1 I 



r ((n - a)/2) 

For the Riesz potential we have the semigroup property 

(2.7) I ai (I a2 f) = I ai + a2 f if ai,a 2 € (0,n) ,s.t. a x + a 2 < n 
and 

(2.8) - A(/ a+2 /) = I a f ifae(0,n-2). 
We need the following simple lemma. 

Lemma 2.2. Suppose that n > 5, a £ [2,n — 3] and -E C ffi" is a set of finite perimeter with 
\E\ = \Bi\. Denote f E — Xb ± — Xe- If the function <fi a : (0, oo) — >• R 



4> a {r):=f {I a f E ){x)<M n - l (x) 

JdB r 

is decreasing, then it holds that 
(i) Q+2 : (0,oo)^K 

0a +2 (r) := / (I a +2f E )(x) dn n -\x) 
JdB r 

is decreasing, 

(ii) 

/ {I a f E ){x)dx<2n{n + 2)-f (I a+2 f E )(x)dx, 

JBx J B\ 

(iii) and 

(I a f E )(x) dy < 2n (I a+2 f E )(0). 



Proof. The proof is standard. Define $ Q : [0,oo) -> 1 as 

:= f (I a f E ){x)dx. 

JB r 

Then it holds 

*a(r)= / / {I a f E ){x)dU n - 1 (x)ds=n\B 1 \ f s n - l 4> a {s)ds. 

JQ JdB B JO 

Since <^> a is decreasing, also i> Q is decreasing. In particular, since lim|, r |_ > . 00 (J Q ,/£)(x) = 0, <fi a and 
<3> Q are non-negative. Therefore, since —A(I a+ 2f E ) = {I a f E ) by (2.8), we have that 

(2.9) ^ +2 (s) = - / (I a+2 f E )(x)dx = ---f (I a f E )(x)dx = --$ a {s) 

nj Bs n JB s n 

and we conclude that </> a + 2 is decreasing and (i) follows. Since lim| a .|_ ) . 00 (J a +2/£;)(x) = 0, (f> a + 2 is 
also non-negative. 

Since <£> Q is decreasing and non-negative we obtain from (2.9) that for every s € (0, 1) it holds 

0' a+2 (s) < ~£*a(l). 

Let r <E (0, 1). Integrate the above inequality over [r, 1] to obtain 

cf> a+2 (l) - d> a+2 (r) = [\' a+2 ( s )ds<-^- f 1 sds = -^-(l-r 2 ). 

Jr J r 

This implies 

(2-10) ^ +2 {r) > a+2 (l) + ^) (1 - r 2 ) > - r 2 ). 

2n 2ra 
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Notice that lining 4>a+2 ( r ) = (Ia+2f E )(0) an d therefore (2.10) implies 

{Ia+2f E ){0) > ^-*a(l) = iri (Ia.f E )w a (x)d X 

In In J Bl 



and (Hi) follows. 

Claim (ii) follows also from (2.10), since 



l 

n-li 



(I a +2f E )(x)dx= / (I a+ 2f E )(x)dU n - x (x)dr 

JO JdB r 

= n|Bi| / r n -Va+a(r)dr 



> ^$ a (l)jf (l-r 2 )r"- x dr 



*a(l). 



2n(n + 2) 



□ 



For the even dimensional cases we also need the fractional Laplacian with a fractional power ^. 
Suppose that / : R" — > R is an integrable function such that 



(i + N 



The ^-Laplacian of / is defined by 



(—A)if(x) := a n [ f & ~J® d£, 

where a n is some normalization constant. The ^-Laplacian can be considered as the inverse of the 
Riesz potential I±, i.e., if / : R 71 — > R satisfies (2.11) we have that 

(2.12) (-A)i(hf) = f. 

The ^-Laplacian can also be obtained by extension problem. Denote R™ + = {(x, z) € M n+1 | 
z > 0}. Suppose that / : R" — > R is a bounded function. Solve the problem 

j -Au(x, z) = for (a;, z) e R^ +1 
| u(a;,0) = /(a;) for x 6 R™. 

Then it holds 

(2.13) -u,(a; J 0) = (-A)3/(x). 

For further details and results about the fractional Laplacian with a general fractional power we 
refer to [2] and [9]. 

Proof of the Theorem 1.1 in dimensions n > 4. The proof begins exactly as in the case n = 
3. We consider the problem in the form (2.1). We may assume that E C R", \E\ = \B\\, is such 
that 

/ „^ f n — 1 , f n — 1 , 
j(E) = / ——dx - / 

Jsx Fl Fl 

We need to show that 

NL(Bi) - NL(E) < P(Bi) - P(E) 

where NL(^) := / R „ f Rn x ^Xf_fa) dydx. Recall that "<" means that the inequality is true up to 
a positive constant. In fact, in the beginning of the proof in n = 3 we have already proven that 

(2.14) NL(Bi) - NL(£) < 2 f v(y) dy, 
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(2-15) v(y)=l <** dx, 



where 

v(y) = f — 

\x - y 
which solves 

-Av = Cn {XB 1 - Xe) 

for some c„ > 0. Define : (0, oo) — > K 

(2.16) 0(r):=/ v(y)dU n - 1 (y). 

JdB r 

For every r > it holds 

(j)'{r) = - f Av(y)dy = — — -f (* Bl - Xs)(y) dy < 
"is,. " is r 

and therefore is decreasing. 

We divide the proof in two cases whether the dimension n is odd or even. 

The case when n > 5 is odd : 

Denote 

Ie = XB 1 - XE- 

Consider the functions /2/B, -?4/b, -?6/b, ■ ■ • , -fn-i/B, where I a is the Riesz potential defined in 
(2.6). In other words 

(XBx ~Xe){x) , , , , _ n-l 



(h k f E ){y) = a n ,2k ! | Bl i dx for k = 1, 2, 



By (2.8) we have that -A(/ 2 fc+2/B) = ^fc/e for every fc = 1, 2, . . . , ^ an d hfE = cr„,2 u , where 
v is defined in (2.15). Define also 02, 04, ■ ■ ■ , <fr n -i '■ (0, oo) — > R such that 

r n—1 
fak(r):=f (hkfEMdW 1 - 1 ^) for fc = l,2,...— -. 

Job,. 2 

Notice that 02 = cr ni 2 0, where is defined in (2.16). Since is decreasing also 02 is decreasing. 
We may use Lemma 2.2 (i) and (ii) iteratively for / 2 /b, J4/, . . . , InsfE to obtain 

(hf E )(y)dy<i (i n -3f E )(y)dy 

Bi J Bi 

and that 0„-3 is decreasing. In other words 

(2.17) [ v(y)dy=[ f <** ~ " *ffi> dxdy < [ f ^'^ dxdy. 

Jb, Jb, Jr- f - 2/| n -/bx Jr" f - yV 

Finally we use Lemma 2.2 (Hi) with a = n — 3, i.e., to functions I n -3f E and to obtain 

5 Bl < y I r^fE)(y)dy < 2n(/ n _i/js)(0). This, in turn, implies 

(XBx-XeX*) ^^ /■ (xbx~Xb)(x) ^ 

F-2/1 Jr™ fI 

The previous inequality and (2.17) imply 

/ (XBx-XB)W & _7(^) 



|ar| n—1 
Combining the above estimate with (2.14) and (2.3), we get 



NL(fli) - NL(E) < 2 f v(y) dy < j(E) < P(E) - P(B X ) 
J Bl 



and the claim follows. 

The case when n > 4 is even: 

Denote 

Ie = XBx - Xb 
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and consider functions I 2 f E , hf E , hfE, ■ ■ ■ , ^1-2 /e, i-e., 

{hklEAV) = <J n .2k / -; ,„_ 2fc cfe for fc = 1, 2, . . . — — 

and functions </>2, 04, . . . , ^> n -2 : (0, oo) — > K such that 

/■ n — 2 

2fc (r) := f (I 2fc / S )(y) cfH"- 1 ^) for fc = 1, 2, . . . — _. 

JdB r 1 

Again <p2 is decreasing and we may use Lemma 2.2 (i) and (ii) iteratively for IiSe^ Ia/e, • ■ • , In-2fE 
to obtain 



(l2f E )(y)dy<f (I n - 2 fE)(y)dy, 
and that </> n -2 is decreasing (and no n- negative) . In other words 
(2-18) f v(y)dy= f f (XBl - X ^f ] dxdy < f (I n - 2 f E )(y)dy. 

J Bi J Bi JR n \ x ~ y\ J Bi 

To conclude the proof we need to show that 



(2-19) / {I n - 2 fEKy)dy<(I n . 1 f E )(0). 

JBx 

The inequalities (2.18) and (2.19) would imply 

v(y) dy < (J n _i/ B )(0) = f iXBl ~* g)(a:) dx < 7 (£) 



and the claim would follow from (2.14) and (2.3). 

To prove (2.19) we use the i-Laplacian. By the semigroup property (2.7) it holds I n _i(f E ) = 
Ii(I n - 2 fE) and therefore by (2.12) we have 

(-A)i(I n _ 1 f E ) = I n _ 2 f E . 
Suppose that u : -> K is the harmonic extension of I n -ifE, i.e., 

j -Au(y, z) = for (y, z) e K™ +1 

\ «(y,0)=J n _i/ B (v) foryeM". 

Then -u z (y,0) = (-A)i(I n -if E )(y) = (/„- 2 /s)(y) for y E f". 

Extend it to M' l+1 by reflection, i.e., define u : M ,l+1 -> K such that 

{■u(y, 2;) = u(y, z) for z > 
u(y, z) = u(y, —z) for z < 0. 

Denote V = (y, z) g R ,l+1 and denote by B r the (n + l)-dimcnsional ball with radius r and 
\Bi\ = uJn+i- Define further fa : (0, oo) — > K such that 

0«(r):=/ fi(y)dH n (y). 

By the divergence theorem we get 

= / (Vf t (r),-)rfH"(r) 

= F - TTT » ( f (Vu(Y),v) dU n (Y) f -u z (Y) dU n (Y)) 

(n + l)uj n+1 r n \J d (B r nRl +1 ) JB r ndRl +1 I 

(2.20) / \ 

Au(Y)dY- f (—A)i(I„-xf E )(y) dy) 



(n + l)uj n+1 r' n 



= -2 7— tv f (I n -2f E )(y)dy, 

(n + l)cj„ +1 
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where ui n is the measure of the n-dimensional unit ball. Since 4> n -2( r ) — fg B (-^n— 2/e)(j/) dH n 1 (y) 
is decreasing and non-negative, also r i— >■ j (I n -2fE)(y) dy is decreasing and non-negative. There- 
fore (2.20) yields 

^(r)<-C n -f {I n - 2 f E )(y)dy, 

J Bi 

for C n = ( n+ i)^ ~ ■ Integrate the above inequality over [R, 1] to obtain 

0u(l) - MR) < ~Cn (j g (I n -2f E )(y) dy^j (1 - R). 

Since limfl'_ 5 .o <p u {R) = u(0) = (7„_i/e)(0) we obtain from the above estimate that 
(1„-i/b)(0) > <fe(l) + C„ / (/„_ 2 / B )(») > (5„ / (I n - 2 f E )(y) dy, 

J B± J B\ 

which proves (2.19). □ 
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